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1. Introduction 

Photodissociation is the dominant removal process of molecules in any region exposed to intense ultra- 
violet (UV) radiation. This includes diffuse and translucent interstellar clouds, dense photon-dominated 
regions, high velocity shocks, the surface layers of protoplanetary disks, and cometary and exoplanetary 
atmospheres. In the simplest case, a molecule ABC absorbs a UV photon into an excited electronic state 
which subsequently dissociates to AB + C. In reality, the situation is much more complex because there 
are many electronic states that can be excited with only a fraction of the absorptions leading to dissocia- 
tion. Also, there are many possible photodissociation products, each of which can be produced in different 
electronic, vibrational and rotational states depending on wavelength. The rate of photodissociation de- 
pends not only on the cross sections for all of these processes but also on the intensity and shape of the 
radiation field at each position into the cloud. Thus, an acccurate determination of the photodissociation 
rate of even a simple molecule like water involves many detailed considerations ranging from its electronic 
structure to its dissociation dynamics and the specifics of the radiation field that the molecule is exposed 
to. 

In this chapter, each of these steps in determining photodissociation rates is systematically discussed. 
§2 reviews the basic processes through which small and large molecules can dissociate following UV 
absorption. Techniques for determining absorption and photodissociation cross sections through theo- 
retical calculations and laboratory experiments are summarized in §3. This section also mentions new 
■ experimental developments to determine branching ratios to different products. Different interstellar and 

| circumstellar radiation fields are summarized in §4, including their attenuation due to dust and self- 

shielding deeper inside the cloud. The formulae for computing photodissociation rates are contained in 
§5, together with references to recent compilations of cross sections and rates. Recent developments in our 
understanding of the photodissociation of CO and its isotopologues are presented in §6. The chapter ends 



6 

CZ5 



with a brief review of our knowledge of the ph otostability of polyc yclic aromatic hydrocarb o ns (P AHs 
Much of this chapter follows earlier reviews by van Dishoeckl |l988l ) and van Dishoeck et al. ( 20061 ) . 



2. Photodissociation processes 

2.1. Small molecules 



The p roces ses by which ph otodissociation of simple molecules can occur have been outlined by Kirbv fc van Dishoeck 



( 19881 ) and Schinke ( 19931 ). A summary is presented in Fig. Q] for the case of diatomic molecules. Similar 



processes can occur for small polyatomic molecules, especially if the potential surface is dissociative along 
one of the coordinates in the multidimensional space. 

The simplest process is direct photodissociation, in which a molecule absorbs a photon into an excited 
electronic state that is repulsive with respect to the nuclear coordinate. Since spontaneous emission back 
to the ground state is a comparatively slow process (typical Einstein-A coefficients of 10 9 s _1 compared 
with dissociation times of 10 13 s _1 ), virtually all of the absorptions lead to dissociation of the molecule. 
The resulting photodissociation cross section is continuous as a function of photon energy and peaks close 
to the vertical excitation energy according to the Franck-Condon principle. The width is determined by 
the steepness of the repulsive potential: the steeper the curve, the broader the cross section. Its shape 
reflects that of the vibrational wavefunction of the ground state out of which the absorption occurs. Note 
that the cross section is usually very small at the threshold energy corresponding to the dissociation 
energy D e of the molecule; thus, taking D e as a proxy for the wavelength range at which dissociation 
occurs gives incorrect results. 

In the process of predissociation, the initial absorption occurs into a bound excited electronic state, 
which subsequently interacts non-radiatively with a nearby repulsive electronic state. An example of such 
a type of interaction is spin-orbit coupling between states of different spin multiplicity. Another example 
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Figure 1. Photodissociation processes of diatomic molecules and corresponding cross sections. From top to 
botto m: direct photodissociation, predissociation and spontaneous radiative dissociation (based on Ivan Dishoeckl 
QUI). 

is the non-adiabatic coupling between two states of the same symmetry. The strength of the interaction 
depends sensitively on the type of coupling and on the energy levels involved, but predissociation rates 
are typically comparable to or larger than the rates for spontaneous emission. The effective photodis- 
sociation cross section consists in this case of a series of discrete peaks, which reflect the product of 
the oscillator strength of the initial absorption and the dissociation efficiency of the level involved. The 
width is controlled by the sum of the radiative and predissociation rates and is generally large (10 12 s _1 , 
corresponding to 15 km s _1 or more in velocity units). 

If the excited bound states are not predissociated, spontaneous radiative dissociation can still be effec- 
tive through emission of photons into the continuum of a lower-lying repulsive state or the vibrational 
continuum of the ground electronic state. The efficiency of this process is determined by the competition 
with spontaneous emission into lower-lying bound states. The photodissociation cross section again con- 
sists of a scries of discrete peaks, but the peaks are not broadened and have widths determined by the 
total radiative lifetime (typically <0.1 km s _1 ). 

Because a molecule has many excited electronic states that can be populated by the ambient radiation 
field, in general all of these processes will occur. As an example, Fig. [2] shows the potential energy 
curves of the lowest 16 electronic states of the OH molecule: predissociation occurs through the lowest 
excited A 2 £ + electronic state, whereas direction dissociation can take place through the 1 2 S _ , 1 2 A 
and coupled 2 and 3 2 H electronic states. However, usually only one or two of these processes dominate 
the photodissociation of a molecule in interstellar clouds. In the case of OH, these are the direct 1 2 S~ 
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Figure 2. Potential energy curves of the OH molecule (reproduced from Ivan Dishoeck fc Dalgarnoiri984f ). 

and 1 2 A channels. However, predissociation through the lower-lying A state is important in cometary 
atmospheres and disks around cool stars, where the radiation field lacks high-energy photons. 

The photodissociation of many simple hydrides like H^, OH, H2O, CH, CH + , and NH proceeds mostly 
through the direct process. On the other hand, the photodissociation of CO is controlled by predissociation 
processes, whereas that of H 2 occurs exclusively by spontaneous radiative dissociation for photon energies 
below 13.6 eV. Whether the photodissociation is dominated by continuous or line processes has important 
consequences for the radiative transfer through the cloud (see §4.6). 

2.2. Large molecules 

For large molecules, the same processes as illustrated in Fig. [T] can in principle occur but they become 
less and less likely as the size of the molecule increases. This is because the number of vibrational degrees 
of freedom of a large non-linear molecule of 3./V — 6 increases rapidly with the number of atoms N. Each 
of these modes has many associated vibrational levels, with quantum numbers v increasing up to the 
dissociation energy of the potential. For sufficiently large N, the density of vibrational levels becomes so 
large that they form a quasi continuum with which the excited states can couple non-radiatively (Fig. [3]) . 
Through this so-called process of internal conversion, the molecule ends up in a highly-excited vibrational 
level of a lower-lying electronic state. The probability is small that the molecule finds its way across the 
multi-dimensional surface to a specific dissociative mode, so the most likely outcome is that the molecule 
ends up in an excited bound level and subsequently relaxes by emission of infrared photons. Alternatively, 
the excited molecule can fluoresce down to the ground state in a dipole-allowed electronic transition, or 
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phosphorescence 



Figure 3. Radiative and radiationless decay of large molecules, illustrating the processes of internal conversion 
(IC), intersystem crossing (ISC), fluorescence, phosphorescence and infrared emission. This example is for a large 
ion such as PAH + with a ground state with one unpaired electron (doublet spin symmetry) denoted as Do, and 
doublet and quartet excited electronic states. Various vibrational levels within the electronic states are indicated. 



it can undergo so-called intersystem crossing to an electronic state with a different spin multiplicity from 
which it can phosphoresce down in a spin-forbidden transition. 

Statistical arguments using (modified) Rice-Ramsperger-Kassel-Marcus (RRKM) theory suggest that 
iV-atom molecules with N > 25 are stable with respect to photod issociation, although this num ber depends 
on the structure and types of modes of the molecule involved ( Chowdarv fc Gruebeld 2008 ). Thus, long 
carbon chains may have a different photostability than ring molecules. PAHs are a specific class of 
large molecules highly relevant for astrochemistry, and are discussed further in §7. Very few gas-phase 
e xperiments exist on s uch large molecules to test these theories. 

Ashfold et al. I (|2010l) have recently pointed out that many intermediate-sized molecules with N rs 10, 



including complex organic molecules such as alcohols, ethers, phenols, amines and N-containing hetero- 
cycles, have dissociative excited states that greatly resemble those found in smaller molecules like H2O 
and NH3. Experiments show that H-loss definitely occurs in these larger systems. 

The ionization potentials of large molecules lie typically around 7-8 eV. Most of the absorptions of 
UV photons with larger energies are expected to lead to ionization but a small fraction can also result 
in dissociation through highly-excited neutral states lying above the ionization threshold or through low- 
lying dissociative states of the ion itself. These branching ratios are not well determined experimentally. 



3. Photodissociation cross sections 

Quantitative information on photodissociation cross sections comes from theory and experiments. Theory 
is well suited for small molecules, in particular radicals and ions that are not readily produced in the 
gas phase. Laboratory measurements of absorption spectra over a wide photon range exist primarily for 
chemically stable molecules, both small and large. In the following, we first review the theoretical methods 
to compute potential energy surfaces and cross sections, and then discuss the experimental data. 

3.1. Theory 

All information about a molecule is contained in the wave function 4', which is the solution of the time- 
independent Schrodinger equation 

HV(x,R) = EV(x,R) (1) 

Here x stands collectively for the spatial and spin coordinates of the n electrons in the molecule and R 
denotes the positions of all N nuclei in the molecule. The total Hamiltonian H consists of the sum of 
the kinetic energy operators of the nuclei a and the electrons i and of their potential energies due to 
mutual interactions. Equation (1) is a (3n + 3A^)-dimensional second order partial differential equation 
that cannot be readily solved, even for small molecules. 
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3.1.1. Born-Oppenheimer approximation 

Because the mass of the nuclei is much larger than that of the electrons, the nuclei move slowly compared 
with the electrons. Most molecular dynamics studies therefore assume the Born-Oppenheimer approxi- 
mation for separating the nuclear and electronic coordinates: 

*(x,B) = * e '(x;R)*" MC (x,R) (2) 

where the electronic energies E el (R) (also called the potential energy curves or surfaces) are determined 
by solving the electronic eigenvalue equation with the nuclei held fixed: 

-H e '* e '(x;R) = E el (R)* ei (x; R) (3) 

Note that now only depends parametrically on the nuclear positions R. Substituting Equation (2) 
into (1) and using (3) gives 

[- E(^ M «) V « + Eel W - E]^ nuc {R) = (4) 

a 

where the sum is over nuclei a with mass M a . In this equation, the assumption V 2 \E ,e '\E ,rmc = ^ el \7^ nuc , 
inherent in the Born-Oppenheimer approximation, has been used. Coupling terms involving V Q ^E' e ' are 
called non-adiabatic terms and their neglect is usually justified. Atomic units (a.u., not to be confused 
with astronomical units or arbitrary units) have been adopted, which have K = m e = e — 1. The unit of 
distance is then 1 a.u. = 0.52918 A (also called Bohr or a ) and the unit of energy is 1 a.u. = 27.21 cV 
(or Hartree). 

Consider for simplicity a diatomic molecule with internuclear distance R. The probability of an elec- 
tronic transition from state i to state / is governed by the magnitude of the electronic transition dipole 
moment D(R), which can be computed from 

D(R) =<*f(r;R)\d\9f(r;R) > (5) 

where the integration is performed over the electron coordinate space and d is the electric dipole moment 
operator in atomic units. 

The photodissociation cross section (in cm 2 ) following absorption from a bound vibrational level v" of 
the ground electronic state into the vibrational continuum k' of an upper state at a transition energy AE 
is then given by 

a v „(AE) = 2.69 x 1(T 18 5 A£| < ^ c (R)\D(R)\^Jf c (R) > | 2 (6) 

where the integration is over the nuclear coordinate R. Here g is a degeneracy factor (equal to 2 for 
n <— E transitions, 1 otherwise) and all quantities are in atomic units. Similarly, the absorption oscillator 
strength between two bound states is 

/„,„„ = \gAE v , v „\ < n? c (R)\D(R)W(R) > | 2 (7) 

Theoretical calculations of photodissociation cross sections thus consist of two steps: (1) calculation of 
the electronic potential curve or surfaces and (2) solution of the nuclear motion under influence of the 
potentials. 

3.1.2. Electronic energies: method of configuration interaction 

Step 1 is the calculation of the electronic potential energy surfaces E(R) and transition dipole moments 
D(R) connecting the excited states with the ground state as functions of the nuclear coordinate R. It 
is important to realize that most aspects of chemistry deal with small energy differences between large 
numbers. For example, the binding energy of a molecule or the excitation energy to the first excited state 
is typically only 0.5% of the total energy of the molecule. Thus, the difficulty in quantum chemistry is 
not only in dealing with a 3n + BN many-body problem, but also in reaching sufficient accuracy in the 
results. For example, the Hartree-Fock method, in which all electrons are treated as independent particles 
by expanding the molecular wave function in Equation (2) in a product of n one-electron wave functions, 
fails because it neglects the 'correlation energy' between the electrons: when two electrons come close 
together, they repel each other. This correlation energy is again of order of a fraction of a % of the total 
energy. 

Over the last decades, different quantum chemical techniques have been developed that treat these 
correlations and there are several publicly available programs. Most of the standard packages, however, 
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such as the popular G AUSSIAN package or packages based on density functional theory (|te Velde et al 



2001 ; Shao et al. 2006), are only suitable and well tested for the ground electronic state and for closed- 



shell electronic structures. More sophisticated techniques based on the CASSCF (complete active space 
self-consistent field) or coupled cluster methods work well for the lowest few excited electronic states 
of a molecule, but these states often still lie below the dissociation energy and do not contribute to 
photodissociation. Accurate calculation of the higher- lying potentials through whi ch photodissociation can 
proceed requires multi-ref erence configuration interaction (MR-CI) techniques (iBuenker fc Peverimhofi 



1974lBuenker et al.lll978l ). for which only a few packages exist (e.g., MOLPRO. IWerner et al-lbOlOtT 



In the CI method, the wavefunction is expanded into an orthonormal set of symmetry-restricted con- 
figuration state functions (CSFs). 

M 

$ e, (x,..x n )=^c s $ s ( X ,..xJ (8) 
l 

The CSFs or 'configurations' generally are linear combinations of Slater determinants, each combination 
having the symmetry and multiplicity of the state under consideration. The Slater determinants are 
constructed from an orthonormal set of one-electron molecular orbitals (MOs; obtained from an SCF 
solution to the Hartrcc-Fock equations), which in turn are expanded in an elementary set of atomic 
orbitals (AOs; called 'the basis set') centered on the atomic nuclei. The larger M, the more accurate the 
results. 

In the multi-reference technique, a set of configurations is chosen to provide the reference space. For 
example, the X 2 II ground state of OH has as its main configuration la 2 2a 2 3a 2 Itt 3 , where a and tt are the 
molecular orbitals. The reference space could be chosen to consist of all possible configurations of the same 
overall symmetry that have a coefficient c s greater than some threshold in the final CI wave function. For 
the OH example, this includes configurations such as la 2 2a 2 3a 1 Aa 1 \iT 3 in which one electron is excited 
from the highest occupied MO (HOMO) to the lowest unoccupied MO (LUMO). Alternatively, an 'active 
space' of orbitals can be designated, e.g. the 2a — 4a and Itt — 2tt MOs in the case of OH, within which 
all configurations of a particular symmetry are considered. The CI method then generates all single and 
double (and sometimes even higher) excitations with respect to this set of reference configurations, several 
hundreds of thousands of configurations in total, and diagonalizes the corresponding matrix. The resulting 
eigenvalues are the different electronic states of the molecule, with the lowest energy root corresponding 
to the ground state. 

The quality of an MR-CI calculation is ultimately determined by the choice of the atomic orbital basis 
set, the choice of the reference space of configurations, and the number of configurations included in the 
final CI. The basis set needs to be at least of 'triple-^' quality (i.e., each occupied atomic orbital Is, 
2p, ... is represented by three functions with optimized exponents £). Also, polarization and Rydberg 
functions (i.e., functions which have higher quantum numbers than the occupied orbitals, e.g. 2p for H, 
3d and 4s for C) need to be added. The number of AOs chosen in the basis set determines the number of 
MOs, which in turn determines the number of configurations. For a typical high-quality AO basis set, the 
latter number is so large, of order 10 8 , that some selection of configurations needs to be made in order 
to handle them with a computer. For example, all configurations which lower the energy by more than 
a threshold value (in an eigenvalue equation with the reference set) can be chosen to be included in the 
final CI matrix. 

Once the electronic wave functions for the ground and excited electronic states have been obtained, 
the expectation values of other operators, such as the electric dipole moment or the spin-orbit coupling, 
can readily be calculated. 



3.1.3. Nuclear dynamics: oscillator strengths and cross sections 

Step 2 in calculating photodissociation cross sections consists of solving Equation (4) to determine the 
nuclear wave functions using the electronic potential curves from Step 1. Often the nuclear wave function 
is first separated into a radial and an angular part. If the angular part is treated as a rigid rotor, the 
radial part can be solved exactly by one-dimensional numerical integration for a diatomic molecule so 
that no further uncertainties are introduced. Once the wave functions for ground and excited states are 
obtained, the cross sections can be computed according to Equation (6). 

For indirect photodissociation processes, the oscillator strengths into the discrete upper levels are 
computed according to Equation (7) for each vibrational level v' . If the coupling between the upper bound 
level with the final dissociative continuum (as computed in Step 1) is weak, first order perturbation theory 
can be used to calculate the predissociation rates k pr in s -1 . The predissociation probability r\ u is then 
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Figure 4. Absorption spectrum of gaseo us H2O illustrating the different electronic states through which pho- 
todissociation can occur (reproduced from iFillion et al.ll2001l ). 



obtained by comparing k pr with the inverse radiative lifetime of the molecule A rad : r/ u — k pr / (k pr + A rad ) . 
In our example of OH, this approximation works well for the calculation of the predissociation of the 
A 2 £ + state for v' ^ 2. If the coupling is strong, as is the case for the OH 2 and 3 2 n states, the coupled 
equations for the excited states (i.e., going beyond the Born-Oppenheimer approximation) have to be 
solved in order to compute the cross sections. 

For (light) triatomic molecules, the calculation of the full 3D potential surfaces and the su bsequent dy 



namics on those surfaces is still feasible , albeit with significantly more efforts involved (e.g.. iKroes et al 



1993t ISchinkd 11993 : Bearda et al. 1995 ). For multi-dimensional potential surfaces, often time-dependent 



wavepacket propagation methods are preferred to solve the dynamics rather than the time-independent 
approach. Besides accurate overall cross sections, such studies give detailed insight into the photodissoci- 
ation dynamics. For example, the structure seen in the photodissociation of H2O through the hrst excited 
A state is found to be due to the symmetric stretch of the excited molecule just prior to dissociation. 

For larger polyatomic molecules, such fully flexible calculations are no longer possible and one or more 
nuclear coordinates need to be frozen. In practice, often the entire nuclear dynamics calculations are 
skipped and only electronic energies and oscillator strengths at the ground state equilibrium geometry 
are computed. Assuming that all absorptions into excited states with vertical excitation energies larger 
than D e lead to dissoc iation (i.e., r\ u = 1), this method p rovides a rough estimate (upper limit) on the 



dissociation rate (e.g., Ivan Hemert &: van Dishoeckl 120081) . This method works because of the Franck- 



Condon principle, which states that the highest cross sections occur at the vertical excitation energies 
(Fig. QJ. Because only a single nuclear geometry has to be considered, the amount of work is orders of 
magnitude less than that of a full electronic and nuclear calculation. 

The overall accuracy of the cross sections and oscillator strengths is typically 20-30% for small molecules 
in which the number of active electrons is at most ^30. Of order 5-10 electronic states per molecular 
symmetry can be computed with reasonable accuracy. Transition energies to the lowest excited states 
are accurate to ~0.1 eV, those to higher states to 0.2-0.3 eV. Thus, quantum chemistry cannot predict 
whether a molecule has an electronic state whose excitation energy matches exactly that of, say, Lyman 
a at 10.2 eV. However, if it has an electronic state with a vertical energy close to 10 eV and absorption 
into the state is continuous, theory can firmly state that there is a significant cross section at Lyman a. 
An example is the 1 2 A state of OH: even with ±30 A uncertainty in position, the cross section at 1215.6 
A is found to lie between 1 and 2 x 10~ 18 cm 2 . 

Virtually all molecules have more electronic states below 13.6 eV than can be computed accurately with 
quantum chemistry. Take our OH example: there is an infinite number of Rydberg states converging to 
the ionization potential at 13.0 eV. However, because the intensity of the radiation field usually decreases 
at shorter wavelengths, these higher states do not significantly affect the overall photodissociation rates. 
Their combined effect can often be taken into account through a single state with an oscillator strength 
of 0.1 lying around the ionization potential. 
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3.2. Experiments 

Laboratory measurements of absorption cross sections have been performed for many chemically stable 
species, including astrophysically relevant molecules such as H 2 0, CO2, NH 3 , and CH4. Most of these 
experiments have been performed at rather low spectral resolution, where the individual ro-vibrational 
lines are not resolved. Figure U shows an example of measurements for H2O. A broad absorption continuum 
is observed between 1900 and 1200 A, with discrete features superposed at shorter wavelengths. Note that 
the electronic states responsible for the absorptions at the shortest wavelengths have often not yet been 
identified spectroscopically. 

The absorption of a photon can result in re-emission of another photon, dissociation or ionization of the 
molecule, and most experiments do not distinguish between these processes. If the photon energy is below 
the first ionization potential and if the absorption is continuous, photodissociation is likely the dominant 
process. If not, additional information is needed to infer the dissociation probabilities r\ u . For the H 2 
example, the broad continuum at 1900-1400 A corresponds to absorption into the A X B\ state, which 
is fully dissociative. However, this is not necessarily true for the higher lying discrete absorptions. An 
example is provided by the case of NO for which fluorescence cross sections have been measured directly 
and found to vary significantly from band to band, with r] u significantly less than unity for some bands 
(|Leelll984h . 

Experiments typically quote error bars of order 20% in their cross sections if the absorption is contin- 
uous. For discrete absorptions, high spectral resolution and low pressures are essential to obtain reliable 
cross sections or oscillator strengths, since saturation effects can easily cause orders of magnitude errors. 

Most broad-ban d absorption spectra date from 30-4 yr ago and have been summarized in various 
papers and books ( Hudson 1971 ; Okabd 1978 ; Lee 1984 ) . An electronic compilation is provided through 
the MPI-Mainz UV-VIS spectral atlas. Over the last two decades, emphasis has shifted to experiments 
at single well-defined wavelengths using lasers, in particular at 1930 and 1570 A. The aim of these ex- 
periments is usually to study the dynamics of the photodissociation process at that wavelength, rather 
than obtaining cross sections. Such studies typically target the lower excited electronic states since lasers 
or light sources at w avelengths <1200 A are not commonly available except in specialized laboratories 
( Ubachs et al.lll994h . While these sophisticated experiments have provided much basic insight into pho- 
todissociation dynamics, the restriction to only one or two wavelengths and lack of cross section data 
makes them of limited value for interstellar applications. 

With the recent advent of powerful synchrotron light sources providing continuous radiation over a large 
wavelength range, the pendulum is swinging back to more astrophysically relevant results. For example, 
photodissociatio n cross sections of molecules such as CH3OH are being measured in the X-ray regime up 
to 100 eV (e.g., iPilling et al.l 120071 1. High-resolution data in the FUV regime up to 13.6 eV are eagerly 
awaited. 

3.3. Photodissociation products 

Inclusion of photodissociation reactions in chemical networks requires knowledge not only of the rate of 
removal of species ABC, but also of the branching ratios to the various products, AB + C, A + BC or 
AC + B. Each of these species can be produced with different amounts of electronic, vibrational and/or 
rotational energy, with the absolute and relative amounts depending on wavelength. None of this detail 
is captured in the standard tabulations, which at best give branching ratios between products integrated 
over all wavelengths or at a single wavelength, and neglect any vibrational or rotational excitation. 

A great number of very elegant modern experiments are now available to explore the fragments, based 
on a range of photofragment tran slational spectroscopy and time-resolved photon-electron spectroscopy 
methods (see Ashfold et al. 201dl for references). Product imaging techniques date back more than 20 
years and characterize both the velocities, interna l energies and angular distributions of the photodisso- 
ciation products (e.g., Chandler fc Houstonlll987l ). 

Veloc ity map imaging is a recent example of this technique well suited for probing photodissociation 
(Fig. El lEppink fc Parked Il997t I Ashfold et aLll2006l ). In this method, a molecular beam containing the 
molecule of interest is created by a pulsed supersonic expansion which ensures internal state cooling. After 
passing through a skimmer and a small hole in the repeller plate lens of the ion lens assembly the beam 
is crossed by a dissociation laser beam to form neutral fragments, which are immediately state-selectively 
ionized by a probe laser beam using resonance enhanced multi-photon ionization (REMPI) . The ions retain 
the velocity information of the nascent photofragments. Velocity map imaging uses a special configuration 
of an electrostatic immersion ion lens (combination of repeller, extractor, and ground electrodes) to ensure 
mapping of the ion velocity independent of its position of formation. After acceleration through the ion lens 
the particles are mass-selected by time of flight on reaching the surface of a two-dimensional (2-D) imaging 



Molecular photodissociation 



9 




Laser 



Figu re 5. Experimental set-up for velocity-map imaging of photodissociation products (based on lEppink fc Parked 
11997ft . 



detector, which converts ion impacts to light flashes recorded by a CCD camera. The two dimensional 
image can then be converted to its desired 3-D counterpart using an inverse Abel transformation. The 
lens effectively reduces the 'blurring' of the images. 

Another variant on this technique is Rydberg H-atom time of flight spectroscopy. Here the nascent H 
atoms produced by photodissociation in the ground n = 1 state are excited into the n = 2 state using 
a Lyman a laser at 1216 A and subsequently to a high Rydberg state with n = 30-80 using longer 
wavelength (e.g., 3650 A) radiation. The neutral Rydberg H atoms then reach a micro-channel plate 
where they are field-ionized and detected. 

A completely different set of experiments has been applied to study the branching ratios of large carbon 
chain molecules. Here high velocity collisions combined with an inverse kinematics sche me have shown 



that the most favorable channels are H production from C„H and C3 production from C„ (jChabot et al 
2010). Otherwise the fragmentation behavior is largely statistical. 



An example of the complexity of the situation is provided by the OH and H2O cases. For OH, the 
potential curves in Fig. [2] show that absorption into the 1 2 £~ state results in ground-state 0( 3 P) 
+ H, but absorption into the higher lying 2 A, 2 n and 2 £ + states in excited 0( 1 D) and 0( 1 S). The 
cross sections into these states are such that about equal amoun ts of Q( 3 P) and Q( 1 D) are prod uced 
in the general interstellar radiation field, but only 5% of 0( 1 S) ( van Dishoeck fc Dalgarnol Il983) . An 
astrophysical confirmation of this p rediction is the detection of the red line of atomic O ( 1 D — > 3 P) at 
6300 A from protoplanetary disks ( Storzer fc Hollenbach 1998t lAcke fc van den Ancker 20061 ). On the 
other hand, for radiation fields dominated by lower photon energies such as in comets exposed to solar 
radiation, the main product is ground-state 0( 3 P). 

H2O is one of the few polyatomic molecules for which the product distribution has been well char- 
acterized as a function of energy. Both theory and experiment have shown that absorption at 6-8 eV 
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Figure 6. Comparis on of the general interstellar r adiation field of Draine (1978) (extended for A > 2000 A using 
the representation of Ivan Dishoeck fc Blackl (|1982f) ) with various stellar radiation fields scaled to have the same 
integ rated intensity from 912-2050 A. The scaled NEXTGEN model radiation field of a B9.5 star (jHauschildt et all 
119991 ) is included as well (dash-dotted). The wavelength ra nges where the photodiss ociation of some important 
interstellar molecules occurs are indicated (reproduced from Ivan Dishoeck etHll2006h . 

into the first A electronic state produces OH in the ground electronic 2 H. state with some vibrational 
excitation but little rotational excitation. In contrast, absorption at 9-11 eV into the second B electronic 
state produces OH in highly excited rotational levels of the A 2 S + excited electronic state, but with low 
vibrational excitation. Absorption into even higher excited states results in O + H 2 or O + H + H prod- 
ucts, rather than OH + H. For other polyatomic molecules, experimental information is spotty at best. 
Early experiments were often performed at high pressures, where the observed products cou ld result from 
subsequent chemical reactions in the gas that would not occur at the low densities in space (|Qkabdll978h . 

Some guess of the most likely products can be made o n the basis of produ ct energies (when does a 
product channel open up?) and correlation diagrams (see Ashfold et al. 2010l for examples), but these 
techniques are usually limited to the lowest-lying channels that may not dominate the overall dissociation. 
Often, the astrochemical databases simply assume a statistical distribution over the various products 
averaged over all photon energies. Note that the information on the product distribution o f species like 



CH3OH is relevant not only for gas-phase chemistry but even more so for ice chemistry (jOberg et al 
20091) . 



4. Astrophysical radiation fields 

4.1. General interstellar radiation field 

Determina tions of the a vera ge intensity of t he int erstellar radiation field (ISRF) fall into two main cate- 
gories (see Black ( 19941 ) and van Dishoeck (1994) for critical reviews up to that date). The first method 
is to estimate the number and distribution of hot stars (O, B, A, ... type) in the Galaxy, use a model for 
the dust distribution and properties to determine the absorption and scattering of the stellar radiation, 
and then su m their fluxes to det ermine the ene rgy density at a typical interstellar locati on. This method 
dates back to Habine ( 1968 ) and Draind ( 19781 ) . and has most recently been revisited by Parravano et al 



(2003). The stellar fluxes used in these models are a combination of observed fluxes of early- type stars 
and extrapolations to shorter wavelengths using model atmospheres. 

The s econd method is t o use direct measurements of the UV radiation from the sky at specific wave- 
lengths (|Henrv et al. I ll980l) . Also in this method, model stellar atmospheres are used to provide informa- 



tion at wavelengths not directly observed. In both methods, direct starlight from early B stars is found 
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to dominate the measured interstellar flux. Specifically, most of the flux at a given point in space comes 
from the sum of discrete sources within 500 pc distance from that point. 

The various estimates of the local ISRF in the solar neighb orhood agree remarkably well, within a factor 
of two. The median value found bvjp arravano et al.l (|2003l ) is within 10% of the mean energy density of 
U = 8x 10~ 17 erg cm" 3 A" 1 of |Drainej (1978J), averaged over the 912-2070 A FUV band (A FUV ). The 
strength of the one-dimensional interstellar radiation field is often indicated by a scaling factor Go with 
respect to the flux F in t he Habing field (-Ffuv = c^fuv^fuv) of 1.6 x 10 -3 erg cm~ 2 s _1 , which is a 
factor of 1.7 below that of iDrainel (|l978[) . Thus, the standard Draine field implies Go = 1.7. Care should 
be taken that photodissociation rates used in astrochemical models refer to the same radiation field. 

If a cloud is located close to a young, hot star, the intensity incident on the cloud boundary may be 
enhanced by orders of magnitude compared with the average ISRF. A well known example is the Orion 
Bar, where the intensity is enhanced by a factor of 4.4 x 10 4 with respect to the Draine field. In contrast, 
at high latitudes where few early-type stars are located, the intensity may be factors of 5-10 lower in 
the FUV band than that of the standard ISRF field. Throughout the galactic plane and over timescales 
of a few Gyr, variations in the local energy density by factors of a few are expected due to the birth 
and demise of OB associations within about 30 Myr. Moreover, the ratio of the highest-energy photons 
capable of dissocia ting H2 and CO and t hose in the more general FUV range may vary by a factor of 2 
in time and place ( Parravano et al.ll2003l) . 



4.2. Stellar radiation fields 

The surface layers of disks around young stars form another example of gas in which the chemistry 
is controlled by photodissociation. The illuminating stars range from late B- and early A-type stars 
(the so-called Herbig Ae/Be stars) to K- and M-type stars (the T Tauri stars). For the latter, if the 
stellar atmosphere dominates the flux, orders of magnitude fewer high-energy photons are available to 
dissociate the molecules compared with the ISRF scaled to the same integrated intensity (Fig. [6]). In 
particula r, the number of photons capable of dissociating H2 and CO and ionizing carbon is greatly 
reduced ( Kamp fc van Zadelhonll2001t Ivan Dishoeck et al.l 12006 ). However, if accretion onto the star is 
still taking place, the hot gas in the accreting column p roduces high-energy photons, including Lyman a, 
which can dominate the FUV flux (jBergin et al.ll2003[) . 



4.3. Lyman a radiation 

Fast shocks (velocities >50 km s _1 ) produce intense Lyman a radiation due to collisional excitation 
of atomic hydrogen in th e hot gas or recombination of ionized hydrogen ( Hollenbach Sz McKed 19791 : 
Neufeld fc Dalgarnolll989T) . Another prominent line in shocks is the C III resonance line at 977 A. Fast 
shocks can be caused by supernovae expanding into the interstellar medium, or by fast jets from protostars 
interacting with the surrounding cloud. Accretion shocks onto the young star mentioned above are another 
example, as are exoplanetary atmospheres. Some molecules, in particular H2, CO and N2, cannot be 
dissociated by Lyman a. Other simpl e molecules such as OH, H 2O and HCN can. A recent summary of 
cross sections at Lyman a is given by Ivan Dishoeck et al.l (|2006l ). 



4.4. Cosmic-ray induced photons 

A dilute flux of UV photons can be maintained deep inside dense clouds through the interaction of cosmic 
rays with hydrogen in the following process: 

H or H 2 + CR -4 H+ or H+ + e* (9) 



H, 



H* 



H2 + ■ 
H, + Yiv 



The energetic electrons e* produced by the cosmic rays excite H2 into the B 1 E+ and G 1 Ii u electronic 
states, so that the FUV emission is dominated by Lyman- and Werner-band UV photons. Higher lying 
electronic states contribute at the shortest wavelengths. As a result, the UV spectrum produced inside 
dense clouds resembles strongly th at of a standard hydrogen lamp in the laboratory. Figure [7] shows the 
UV spectrum computed following iGredel et al. I (Il987t Il989h for e*=30 eV and either all H 2 in J=0 or 
distributed over J=l and J=0 in a 3:1 ratio. Note that the precise spectrum does depend on the H2 level 
populations and o/p ratio. Because of the highly structured FUV field and the lack of high-resolution 
cross sections for most species, the resulting photodissociation rates may be quite uncertain. On the other 
hand, the large number of lines has a mitigating effect. 
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Figure 7. Cosmic-ray induced spectrum of H2 assuming that all H2 is in J—0 (left) or distributed over J=l and 
J=0 in a 3:1 ratio (left). The spectrum shows the number of photons per 1 A bin for initial 30 eV electron s and 
are normalised to the total number of ionizations (figure by Gredel, updated from lGredel et al.|[l987l . Il989l ). 



4.5. Dust attenuation 



The depth to which ultraviolet photons can penetrate in the cloud and affect the chemistry depends on 
(i) the amount of photons at the boundary of the cloud; (ii) the scattering properties of the grains as 
functions of wavelength; (iii) the competition between the atoms and molecules with the grains for the 
available photons, which, in turn, depends on the density of the cloud; and (iv) possible self-shielding of 
the molecules. 

The absorption and scattering parameters of the grains determine the penetration depths of photons. 
The ultraviolet extinction curve shows substantial variations from place to place, especially at the shortest 
wavelengths in the 912-1100 A region, which can grea tly affect the CO and CN photodissociation rates 
( Cardelli fc Savage! 19881 : van Dishoeck fc Blacklll989l) . Similar variations are expected to occur for the 
albedo vo and the asymmetry parameter a as functions of wav elength, but little is still known about these 
parameters at the shortest wavelengths ([Roberge et al.lll991l ). 

There is considerable observational evidence that dust grains in the surface layers of protoplanetary 
disks have grown from the typical interstellar size of 0.1 [ivct to at least a few /im in size or more. These 
large grains e xtinguish the UV radiation much less, allowing photodissociation to penetrate deeper into 
the disk (e.g., van Dishoeck et aT1l2006 : Vasvunin et al. 2011 ). 



4.6. Self- shielding 

Because the physical and chemical structure of cloud edges is controlled by the photodissociation of H2 and 
CO and photoionization of C, it is important to pay particular attention to the details for these molecules. 
CO is discussed in §6. The photodissociation of H 2 occurs by absorptions into the Lyman and Werner 
bands at 912-1100 A followed by spontaneous emission back to the vibrational continuum of the ground 
state ( Dalgarno fc Stephens 197Cll . see also Fig. Q]). The absorption lines have typical oscillator strengths 
/ w 0.0001-0.03, and on average about 10% of the absorptions lead to dissociation. The strongest lines 
become optically thick at H2 column densities of 10 14 -10 15 cm~ 2 for a Doppler broadening parameter b 
of 3 km s , implying that molecules lying deeper into the cloud are shielded from dissociating radiation 
because all relevant photons have been absorbed at the cloud edge ('self-shielding'). The photodissociation 
rate of H2 is about 4.5 x 10~ n s _1 in the unattenuated interstellar radiation field, corresponding to a 
lifetime of about 700 yr. Inside the cloud, the lifetime becomes longer by several orders of magnitude 
because of the self-shielding process. The remaing 90% of the absorptions are followed by UV fluorescence 
back into the bound vibrational levels of the H 2 ground state. 

Photoionization of atomic carbon has a continuous cross section of about 10~ 17 cm 2 over t he 912- 
1100 A region. Thus, for column densities N(C) > 10 17 cm 2 , carbon starts to self-shield (jWernerl 
1970). Moreover, the saturated absorption bands of H and H2 over the same wavelength range remove a 



considerable fraction of the ionizing radiation ('mutual shielding'). 
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5. Photodissociation rates 

The photodissociation rate k p d of a molecule by continuous absorption is given by 



k 



pd 



a(X)I(X)dX bT 1 , 



(10) 



where a is the cross section for photodissociation in cm 2 and I is the mean intensity of the radiation 
in photons cm~ 2 s™ 1 A -1 as a function of wavelength A in A(Fig. [6]). For the indirect processes of 
predissociation and spontaneous radiative dissociation, the rate of dissociation by absorption into a specific 
level of a bound upper state u from lower level I is 



uline 
^pd 



;XlefueVuXeI(Xui) 



(11) 



where f U £ is the oscillator strength, rj u is the dissociation efficiency of the upper level (between and 1), 
and xi is the fractional population in level t. 

The determination of the cross sections and oscillator strengths from theory and experiments has been 
discussed in §3. Cross se ction databases, including references to the original da ta, are available from 
Ivan Dishoeck et al.l (|2006h at www.strw.leidemmiv.nl/^ewine/photo. and from Huebner et al. (1992) 
at amop . space . swri . edu aimed at comctary species. For more complex species of atmospheric interest, a 
large compilation is available through the MPI-Mainz UV-VIS spectral atlas at www. atmosphere .mpg. de/ 
enid/2295 and www. science-softcon.de/spectra. 

Photodissociation rates as functions of dept h into a cloud using the Draine (1 19781) radiatio n field have 



been prese nted for a variety of molec ules by Ivan Dishoeck ( 19881 ) and Roberge et al. ( 1991 ) , and most 



recently by Ivan Dishoeck et ajj (|2006l ). The latter paper and associated website also summarize rates in 
cooler radiation fields such as appropriate for late-type stars. The depth dependence due to continuum 
extinction can be represented by a simple exponential function, with exponents that vary with grain 
scattering properties and that depend parametricall y on cloud th i ckness . A list of cosmic-ray induced 
rates for use in dense cloud models has been given bv lGredel et al. I (Il989l) . 



6. Photodissociation of CO and its isotopologues 

CO is the most commonly observed molecule in interstellar space and used as a tracer of molecular 
gas throughout the universe, from local diffuse clouds to dense gas at the highest redshifts. Additional 
impetus for a good understanding of its photodissociation processes comes from recent interpretations 
of oxygen isotope fractionation in primitive meteorites, which are thought to originate from isotope 
selective photodissocia tion of CO in the upper layers of the disk out of which our solar system formed 
dLvons fc Young|l2005l) . To model these processes, information on the electronic structure of all the CO 
isotopologues, including those with 17 O, is needed. 

CO is an extremely stable molecule with a dissociation energy of 11.09 eV, corresponding to a threshold 
of 1118 A. It took until the late 1980s to establish that no continuum absorption occurs longward of 912 A, 
but that CO photodissociation is domin a ted by line absorpt ions, most of which are strongly pre dissociated 
(jLetzelter et al.lll987l: [stark et al.|[l99lh . IViala et al.l (|l988l ) and lvan Dishoeck fc Blackl (| 19881 ) used these 
early laboratory data to build a detailed model of the CO photodissociation rate into an interstellar 
cloud. The depth dependence of the CO photodissociation is affected not only by self-shielding, but also 
by mutual shielding by H and H2 because these species absorb in the same wavelength region. 12 CO, in 
turn, can shield the less abundant 13 CO, C 18 and C 17 species, with the amount of shielding depending 
on the wavelength shifts in the absorption spectrum. Thus, a complete numerical simulation of the entire 
spectrum of 12 CO, its isoto pologues, H2 and H is required to co rrectly compute the attenuation at each 
depth into a cloud (see also lLee et al.lll996t iLe Petit et al.ll2006t ). 

In the 20 years since these first models, a steady stream of new laboratory data has become available 
on this key process. In the late 1980s, many of the bands recorded in the laboratory were for 12 CO only 
and had no electronic or vibrational designations, so that simulation of the isotopologue spectra involved 
many guesses. In particular, the magnitude of the isotopologue shifts depends sensitively on whether v' 
in the upper state is zero or not. Moreover, predissociation rates were generally not known and were 
simply assumed to be unity. Hi gh- resolution spectra for many important bands of * 2 CO, 13 CO, C 18 
and 13 C 18 were measured by (jEidelsberg et al.lll992[ |2004 l200a IStark et aT1ll99ll) with synchrotron 
sources. The fact that the experimental oscillator strengths have now bee n reconciled with v alues inferred 

H2003h. 



from astronomical measurements leaves little doubt on their accuracy dSheffer et al 



Ultra-high- 



resolution spectra of selected states have been obtained with VUV lasers (jUbachs et al.lll994j ). These 
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Figure 8. Photodissociation and photoionization rates of PAHs as functions of number of carbon atoms for the 
general interstellar radiation field. Solid: loss of C, dotted: loss of H; dash-dotted: i oniza tion. Note the stability 
of the larger PAHs against dissociation. Figure made using data from IVisser et all <|2007l ) following |AU ain et al.l 
111996ft . 

exquisite data provide not only positions down to 0.003 cm" 1 accuracy, but also measure directly the line 
widths and thus the predissociation probabilities. This is especially important for the C X S + and E 1 n 
states responsible for mo st of the isotope-selective effects ( Cacciani et al. 19951 120011: lUbachs et al.lbood 
ICacciani fc Ubachsll2004l) . 

IVisser et al] (]2009f ) used the recent laboratory data to develop a new model for the CO isotope-selective 
photodissociation including the 17 isotopologues, and applied it to interstellar clouds and disks around 
young stars. They also computed shielding functions for a much larger range of CO excitation tempera- 
tures and Doppler broad enings. Although the overall rate has changed by only 30%, from 2.0 to 2.6 x 10~ 10 
s _1 in the iDrainel (1978) field, the depth dependence differs significantly from earlier work. In particular, 
the isotope selectiv e effects are diminis hed by a factor of three or more at temperatures above 100 K. 
Note, however, that IVisser et al. I (l2009h still had to make many assumptions on line positions, oscillator 
strengths and predissociation rates for the minor isotopologues, since high-resolution spectra of these 
minor species have not yet been measured or published. Unexpected differences of up to an order of mag- 
nitude have been found between oscillator strengths of 12 CO and 13 CO for the same bands, complicating 
extrapolations from the main isotopologues. These differences likely arise from mixing between various 
electronic states that depend sensitively on the relative location of the energy levels and may thus differ 
for a specific ro-vibronic v', J' level of a particular isotopologue (see lStark et al.ll2008l for the case of N2) 



Ano ther set of experiments on the CO isotopologue photodissociation has been carried out bv lChakraborty et al 
( 2008h using the Berkeley synchrotron source. In contrast with the Paris and Amsterdam results, these 
data are at comparatively low spectral resolution and were analyzed much more indirectly following a 
set of subsequent chemical reactions. The data were interpreted to imply different predissociation prob- 
abilities of individual lines of isotopologues due to near-resonance accidental predissociation processes. 
In this interpretation, the isotope-selective effects found in meteorites would not require self- or mu- 
tual shielding processes. While accidental pre-dissociation may indeed occur, these experiments and their 
conclusions have been challenged on many grounds and are not supported by the higher r esolution data 
cited above (|Federman fc Yountdl^OOl iLvons et aLll2009t IyTii et alj|2009l: IVisser et alj(2009h . This episode 



demonstrates that there is no substitute for high-quality molecular physics data in which individual un- 
saturated lines are recorded, in order to draw astrophysical conclusions. 



7. Photostability of PAHs 

PAH molecules are a class of large molecules that are very stable against photodissociat ion in the gen - 
eral interstellar medium, explaining their ubiquitous presence in galaxies (see review by iTielensi 2008) . 
However, in clouds exposed to very intense radiation, such as the inner regions of protoplanetary disks 
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Figure 9. Overview of th e PIRENE A exper iment to study the photostability and fragmentation pattern of PAHs 
exposed to UV radiation ( Joblin et al. 2002). 

and the immediate environments of active galactic nuclei, the smaller PAHs can be destroyed within the 
lifetime of the source. As discussed in §2.2, the dissociation process of a large molecule involves a multi- 
step process, in which the molecule first absorbs a photon into an excited electronic state, then decays 
non-radiatively to high vibrational levels of the ground state, and finally finds a path to dissociation. 
The rate of this process thus depends not only on the initial absorption rate, but also on the competition 
with other processes during these steps. In the first step, the absorption of the photon can also lead to 
emission of an electron through ionization or photodetachment with probability r) cm . In the final step, 
there is competition betwee n dissociation and cooling by infrared emission. The PAH photodissociation 
rate thus becomes (following ILeger et al. 19891 : Allain et al. 1996 ; Le Page et al.ll200l ; Habart et al. 2004 ; 
Visser et al1l2007l) : 



/(!-,« 



I 



where r/diss is the yield of a single dissociation process given by 



Vdh 



Here fcdiss.x is the rate for a particular loss ch annel X (e.g., H r emoval), fcdiss is the sum over all such 



channels, and &ir is the infrared emission rate. ILeger etaL identified loss of H, H 2 , C, C2 and C3 

as possible channels, with the rates determined by the RRKM quasi-equilibrium theory according to 

p(E int — E ^x) 



Miss.X 



= A 



x- 



p{E int ) 

where Ax is the pre-exponcntial Arrhenius factor for channel X and p(E) is the density of vibrational 
states at energy E. Dissociation only occurs if the internal energy of the molecul e E-^t exceed s the c ritical 
energy E a for a particular loss channel. Values for E 0l x were summarized bv IVisser et al.l ( 20071 ). The 
above formulation needs to be modified for clouds exposed to very intense radiation fields (Go > 10 4 ), 
when multiphoton events start to become significant, i.e., the PAH molecule absorbs another UV photon 
before it has cooled down completely. This can significantly increase the dissociation rate. 

To compute the dissociation rates, values of <7 a b s as a function of wavelength are essential. The data used 
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in astrophysical mode ls largely come from models developed bv lLi fc Drain! (|200lh . based on experiments 
bv l Joblin et al. ( 1992). T here is also lim ited experimental informatio n on the loss channels, as summarized 
by Le Page et al. ( 2001 ). For example, Jochims et al. ( 1994 [1999^ have measured the photostability of 
small PAHs of various sizes and shapes against H-loss. Fig. [8] compares the H-loss channel quantitatively 
with that of other loss channels and illustrates the competition with ionization. 

In a novel experiment, called PIRE NEA (Fig. Efl), th e spectroscopy and photodissocia tion of PAHs 
into various fragments can be studied ( Joblin et all 2002 . 2011 : Useli-Bacchitta et al. 2010h . Specifically, 
gas-phase PAH cations are produced by laser irradiation of a solid target and then guided into an ion 
cyclotron resonance cell, where they are trapped in a combined magnetic and electric field. The stored ions 
are irradiated by a xenon lamp (2000-8000 A) and the products are analyzed by Fourier-transform mass 
spectrometry. Ions of interest can be selected and isolated by selective ejection of other species. Initial 
re sults show that irradi ation leads predominantly to atomic hydrogen atom loss, consistent with the work 
of I Jochims et al. ( 1994), but los s of C2H2 and perhaps even C4H2 is also seen from dehydrogenated PAHs 
(jUseli-Bacchitta fc Joblm1l2007t ). 



8. Summary 

This chapter has summarized our understanding of basic photodissociation processes and the theoretical 
and experimental approaches to determine cross sections, product branching ratios and rates for astro- 
physically relevant species. The demand for accurate photorates is likely to increase in the coming years as 
new infrared and submillimeter facilities open up the study of molecules in the surfaces of protoplanetary 
disks and exoplanetary atmospheres. Critical evaluation of the photodissociation data for these different 
environments are needed, since rates or cross sections determined for interstellar clouds cannot be simply 
transposed to other regions. 
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